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Tiivistelma The ttN forward scattering data are analyzed using an expansion method, where the invariant 
amplitudes are represented by expansions satisfying the forward dispersion relations. The experimental 
errors of the data are taken into account through the covariance matrix of the coefficients of the expansions 
in a careful error analysis. From the results, some coefficients, c^q, of the subthreshold expansions have 
been calculated with proper error bars. 
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1 Introduction 

Forward dispersion relations are a special case in ttN scat- 
tering analysis, because in forward scattering the optical 
theorem provides a direct connection to total cross sec- 
tion data. An expansion method provides a tool for guar- 
anteeing that the forward dispersion relations are satis- 
fied. Since the last forward analysis with the expansion 
method [1] , there has not been very much experimental ac- 
tivity in the forward ttN scattering. In particular, there are 
no new total cross section data except the very high ener- 
gy measurements of the SELEX collaboration 2J . Howev- 
er, there is new information at the physical threshold from 
the very precise pionic hydrogen experiments [3], and, at 
low energy, from some new integrated cross section mea- 
surements [4l[5] which are filling the gap between the phys- 
ical threshold and the first total cross section data points. 
Both of these have a direct impact on the subthreshold 
expansion. 

The aim of the present article is to construct invariant 
amplitudes at i = to constrain a phase shift anal- 
ysis with fixed-t analyticity. Furthermore, information on 
Im C~ (w) can be used to study the Goldberger-Miyazawa- 
Oehme sum rule (6l[7]- The expansion method is briefly 
explained in sect. [21 the experimental input is described 
in sect. |3l the minimizations are discussed in sect. |4] to- 
gether with the evaluation of the subthreshold expansion 
coefficients. In sect. [S] the conclusions are drawn. 



2 The expansion method 

The isospin even pion-nuclcon C-amplitude satisfies the 
forward dispersion relation [S] 
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Here v is the crossing antisymmetric Mandelstam vari- 
able V = is — u)/'im = (jj + t/Am with to denoting the 
pion total laboratory energy, is the ttN coupling con- 
stant, vb = —n^/2m and C~^{iy = 0,i = 0) is the sub- 
traction constant. Everywhere we denote the proton mass 
by TO and the charged pion mass by fi. In practice, equa- 
tions like ll]) are difhcult to use as constraints, because 
the principal value integration is difhcult to handle with 
experimental input. In particular, integrating over experi- 
mental data is unreliable and propagating the errors of the 
input to the error bars of the output is very difficult. In- 
stead, we express the isospin even and odd amplitudes us- 
ing Pietarinen's versions of the expansions [MTU] at i = 



C±(z.) = C±(z.) + i7±(Z)Vc±[Z( 
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where the base function Z{h') is a conformal mapping 
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which maps the physical cut to the upper half of the unit 
circle in such a way, that the threshold is mapped to (1,0) 
and the infinity to (—1,0). The threshold of the cut in the 
forward direction is vt = [t. and the parameter a controls 
which energy is mapped to (0, 1); there the pion laborato- 
ry momentum piab — a. The numerical value of a in ([3]) 
is not crucial, but it is fixed to a = 0.72 GeV, which 
seems to give the most rapid convergence ^lOj. By using 
Z{i') as the base function, the analyticity structure of the 
invariant amplitudes is a built-in feature and not only a 
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constraint, so the resulting amplitudes will satisfy the dis- 
persion relations exactly. In the expansion method there 
is also the great advantage, that experimental input for 
both real and imaginary parts with their errors can be 
used simultaneously. The pole terms piTT| 
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1 



1 
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are treated separately and the assumed high energy be- 
haviour in the forward direction is taken care of by the 
functions H^[Z{u)) [9] 



^ l + {rilog[r2(l + Z)]f 
^ ' 1 + Z 

H-{Z) = Amv{l + Zf-"'" , 
where ri — 0.144, r2 — 44.154 , and cxr 



(5) 



0.48. 



The number of terms in eq. ((2|), i.e. iV, was taken to be 
40 in the earher work [inifTT] . but here we take N = 100, 
which should fully guarantee that the truncation error is 
negligible. 

An essential ingredient of the expansion technique is 
the convergence test function (CTF) [10T1JT2|. It is not 
sufficient to fix the coefficients in the expansion at 
fixed-t by fitting to data, but, in addition, to guarantee 
the smoothness of the invariant amplitudes an addition- 
al term in the sum is needed. The convergence test 
function part takes the form [10] 
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where the weights are 
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Here and Ac^ denote the expansion coefficients and 
their corresponding errors in the isospin even and odd 
C-amplitudes at the x^-minimum, i.e. at the best fit to 
the data and to the contraints. So, the minimum has to 
be roughly known before the final CTF weights can be 
calculated. This leads to an iterative minimization. 



3 The input 

As input for the fit, we used total cross section data [131 
[14], integrated cross section measurements |4l[5], real-to- 
imaginary ratios [T3lll4j . real parts of the isoscalar D- 
amplitude 15,l61[T7]. the s-wave scattering length a^-p 
from pionic hydrogen experiments [3l[7] and the scatter- 
ing length a^+p from discrete phase shift analysis [7]. To 
begin with, our full forward data base contains 1098 data 
points in 142 data sets covering the laboratory momenta 




Kuva 1. Electromagnetic (solid line) and Z\-splitting correc- 
tions (dashed line) for o"J-* . 




Kuva 2. Electromagnetic correction for crj°p. 



from 0.077 GeV/c to 640 GeV/c in addition to the thresh- 
old values. 

The electromagnetic effects were removed from the 
total cross sections and from the real parts of the D- 
amplitude by the Tromborg method [18] for the labora- 
tory momenta piab < 725 MeV/c. The corrections have 
been published only up to 655 MeV/c, but here we em- 
ploy a smooth extrapolation up to 725 MeV/c. In order 
to apply the Tromborg method, an existing partial wave 
solution is needed. We used the KA84 solution T9], but 
the corrections are practically unchanged, if one chooses 
to take the KH80 [20] or the FA02 solution [21] instead. 
The Z\-splitting was treated by using the P33 phase shift 
differences from ref. [22] , which are very similar to the ear- 
lier P33-corrections of Bugg [231124] . Bugg gives the correc- 
tions for a discrete set of momenta frompiab = 183 MeV/c 
to piab = 408 MeV/c, while with ref. [22] it is possible to 
treat all data up to piab = 725 MeV/c. The corrections 
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cr"""? — ^'^'^■^'^'^'^ applied to the measured to- 
tal cross sections to obtain the isospin invariant hadronic 
ones are displayed in figs. [T] and [2l After applying these 
corrections, the data were assumed to be purely hadron- 
ic and isospin invariant. That is reasonable, because any 
remaining effects are expected to be considerably smaller 
than the experimental errors. Above piab = 725 MeV/c 
the one photon exchange picture with the Coulomb phase 
and form factors was assumed to be applicable. For the pi- 
onic hydrogen results xPT-based electromagnetic correc- 
tions |25j were applied to extract the hadronic quantities. 
To make use of the integrated differential cross section 
data [HE], they were corrected to hadronic cross sections 
integrated over the whole angular range by adding correc- 
tions calculated with the KA84 solution. 

The error bars of the total cross sections of Carter 
et al. [26| were modified by adding the errors due to the 
0.25% uncertainty in the beam momenta, as explained by 
Bugg [27]. Also, the corrections adopted by Giacomelh [28] 
were applied to the total cross sections of Citron et al. . 

The forward data alone are not enough to stabilize 
the low energy behaviour, i.e. the energy range from the 
threshold up to the first total cross section data point. 
In order to stabilize it without introducing any bias from 
earlier solutions, we used our current partial wave solution 
to constrain the momentum range 20 MeV/c < puh < 
155 MeV/c. The details of the partial wave solution will 
be published elsewhere [30] . 



4 Results 

The coefficients of Pietarinen's expansions ^ were 
fixed in a x^-minimization using the program MINUIT [3T1 
I32j . The actual minimization was carried out three times, 
because the convergence test function method depends on 
the previously determined minimum. The data were al- 
lowed to float inside the quoted systematic errors and the 
floating factors were searched simultaneously with the co- 
efficients of Pietarinen's expansions. In the process of the 
minimization we had to discard six data sets, which were 
too discrepant even after the renormalization: the ir^p to- 
tal cross section sets of Devlin et al. [33], two Tr+p total 
cross section sets of Brisson et al. [M1I55] . the 7r+p total 
cross section set of Ignatenko et al. [36] and the Tr+p total 
cross section set of Lindenbaum et al. [37] . In addition to 
these, we had to discard the three lowest data points of 
Davidson et al. [38j . After excluding these data, we were 
left with 967 data points in 136 sets and the average 
per data point was 1.47. The differences between the da- 
ta and the total cross sections calculated from the fit are 
plotted in figs. [3] and |4] for the laboratory momentum range 
Piab =150 - 500 MeV/c. At higher energy the results are 
in good agreement with the earlier work of the Karlsruhe 
group [T]. 

Inside the Mandelstam triangle, it is useful to for- 
mulate the ttN amplitudes in terms of the subthreshold 
expansion with the pseudovector Born terms subtracted. 
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Kuva 3. Difference plot of the total tt~p cross section. The 
dashed lines give the error band of the fit. The Carter data [26] 
are marked with crosses, the Pedroni data [39] with bars and 
the Davidson data [38] with circles. 
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Kuva 4. Difference plot of the total vr+p cross section. For 
the explanation of the symbols for the data, see the caption of 

fig-H 



cients are the coefficients in the expansions based on cross- 
ing M 



C-{,.,t)/,. = J2c-^,y'-f^, 



(8) 



N .pv ' 



The subthreshold expansion coeffi- 



where the isospin odd amplitude is divided by v in order 
to get a crossing even quantity. One gets the subthreshold 
parameters c^q by Taylor expanding Pietarinen's repre- 
sentation ([2]) with the pseudovector Born term subtracted 
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Taulukko 1. The coefficients of the subthreshold expansion of in natural units (powers of ^^^). In the error estimate the 
first part is the statistical error, and the second part is the combination of the uncertainty in the coupling constant and the 
effect of conflicting data sets. They should be added linearly in order to get the total error. The Karlsruhe results are from 
table 2.4.7.1. of ref. [8] (N.B. c^p = d±,). 





Present analysis 


Karlsruhe 


SM99 40j 


c+ 

^00 


-1.201" ±0.004 ±0.03 


-1.46±0.10 


-1.26±0.02 


^10 


1.119±0.001±0.002 


1.12±0.02 


1.11±0.02 


'-20 


0.2015±0.0005±0.0008 


0.200±0.005 


0.20±0.01 


^30 


0.0568±0.0003±0.0001 







The value resulting from a calculation with isospin invariance. 



Taulukko 2. The coefficients of the subthreshold expansion of C~ jv in natural units (powers of The error estimates 

are displayed in the same manner as in table [T] except for Cqq , where the dependence and the effect of confiicting data are 
separated. 





Present analysis 


Karlsruhe 


SM99 go] 




1.41±0.002'' ±0.05'' ±0.01" 


1.53±0.02 


1.43±0.01 




-0.167±0.001±0.001 


-0.167±0.005 


-0.16±0.01 


^20 


-0.0388±0.0004±0.0005 


-0.039±0.002 


-0.04±0.01 


C30 


-0.0092±0.0002±0.0001 







''The statistical error. 

''The dependence, see sect. |4l 

'^Due to confiicting data. 
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The numerical values of the lowest coefficients are given 
in tables [T] and [2] together with the Karlsruhe values [8] 
and the GWU/VPI SM99 values 

The statistical errors of the subthreshold expansion 
coefficients ([9|) are calculated by the standard way, i.e. 



EdCnO dCnO 
—Vki, 



k.l 



dck dci 



(10) 



where Vki is the covariance matrix of the Pietarinen coef- 
ficients Cfc, calculated by MINUIT. In Pietarinen's method, 
the minimized function is not a pure ^^-distribution, but 
the combination of a ^^-distribution and the convergence 
test function, which generally makes the minimum steep- 
er and gives too optimistic error bars. In practice 9.5% 
of the probability distribution at the minimum is due to 
CTF. Therefore, in order to get the proper error bars, we 
increased the statistical errors of eq. (fTO]) by 10%. 

Two sources of systematic errors are studied explicitly. 
The effect of the uncertainty in the coupling constant was 
estimated by making calculations with various valuesEI in 
the range p = 0.075 ± 0.002 7J. The conflicting data sets 
are causing another systematic effect, which was estimated 
by making the analysis with different subsets of the data. 



^ The pseudoscalar coupling constant g is related to the pseu- 
dovector coupling f hy g^ = 4tt f'^ {2m/ fi)^ . 



The combinations of these effects are displayed in tables 
[T]and[l 



The coefficient Cqq can be written as 



1 



4 — fj, 



(11) 



where a; = ^i/m, a^^ is the isoscalar s-wave scattering 
length and is the integral 



V 

TT 



cr+(fc)dfc 



(12) 



Here a+ is the average a+ = {a^°^ ± a'^f )/2. The un- 

certainty in Cqq is mainly due to the first term of eq. (Ilip . 
The value of the integral J+ is quite stable, we obtain 
J+ = 1.459 ±0.005 /Lt~^, if the integrations are performed 
as in ref. . The result can be compared with the Karls- 
ruhe value [H] J+ = 1.478 ± 0.010 yi-^. If the values 
for a^rip derived in ref. [7] are used and isospin invari- 
ance is assumed, we obtain Oq, = 0.0085 ± 0.0016 11^^ 
which gives the value in table [T] However, employing the 
isovector scattering lengths a^j^ derived from the width 
measurements of pionic hydrogen to fix the isoscalar scat- 
tering length from a^_^ = a.,^-p — a^j^ would give a range 
of values for cJq extending from —1.17 to —1.41 . 
The detailed value depends on which data one wants to 
fit, e.g., the Denz et al. Re_D^ data |17j favour a more 
negative cJq. 

To fix Oqj^ it is also possible to make use of the mea- 
sured level shift in pionic deuterium [42j . For recent eval- 
uations of the 7r"(i level shift, see e.g. ref. [43ll44p 45 p 46j . 
MeiBner et al. [13] include isospin breaking corrections to 
0{p^) in chiral perturbation theory and give the result 
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0.0015 ± 0.0022 which would lead to c, 



— 1.30 /i ^, a result well within the above range of values. 

In table [21 the explicit dependence of the param- 
eter C(^Q is given separately as well as the effects due to 
conflicting data sets. For all other parameter values the 
coupling constant dependence, the conflicting data effects 
and the statistical errors are added linearly. Beyond cJq 
all the other parameters are very stable. 

If one uses the value of the coupling constant from 
Ericson et al. [47ll48]. g'^/^ir = 14.07 ± 0.17 instead of the 
value of ref. [7], the value of the first coefficient of the 
isospin odd expansion will be Cqq — 1.49 ± 0.03 . The 
changes in the other coefficients are small. 



5 Conclusions 

The resulting expansions give smooth forward isospin even 
and odd C-amplitudes, which can be used as a starting 
point for a phase shift analysis with fixed-i constraints as 
well as constraints for a discrete phase shift analysis [149 . 
Also, when the Z\-splitting corrections have been taken 
into account, the hadronic Tr^p total cross sections can be 
calculated and used to evaluate the Goldberger-Miyazawa- 
Oehme sum rule i7|. 

The resulting forward subthreshold expansion coeffi- 
cients are in excellent agreement with the earlier evalua- 
tions of the Karlsruhe group except for c^q, where the ef- 
fects due to the coupling constant and the threshold am- 
plitudes, obtained from pionic hydrogen experiments, are 
important. 
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